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We describe how complications encountered in the extraction of V us from flavor-breaking hadronic 
t decay sum rules, in particular those associated with the slower than previously anticipated conver- 
gence of the relevant D = 2 OPE series, can be successfully dealt with through judicious choices of 
sum rule weight. Problems with the previously proposed "(0, 0) spectral weight" version of the anal- 
ysis, and the much improved situation produced by the use of a number of alternate, non-spectral- 
weights, are demonstrated. The method is shown to promise a sub-O.OOlO-accuracy determination 
of Vus once analysis of the improved strange spectral data from BABAR and BELLE is completed. 
Connections to the r versus electroproduction puzzle in the evaluation of the leading order hadronic 
contribution to (g — 2) M are also discussed. 
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where y T = s/m 2 ., Vij is the flavor ij CKM matrix el- 
ement, and Sew is a short-distance electroweak correc- 



tion. Analogous spectral integrals, R 



V/A;ij ' 



are obtained 



from Rv/A;ij by rescaling the integrand with (1 — y T ) k y r ^ 
before integration. In what follows we denote by R^(so) 
the weighted spectral integral from threshold to s = sq < 
m 2 involving the analytic weight, w(s), and either of the 



correlators IT^^^s) or sll^y^.^s). All spectral inte- 
grals discussed above satisfy the general finite energy sum 
rule (FESR) relation 
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and hence have an alternate, OPE representation for suf- 
ficiently large sq. The FESR's corresponding to Rv/A\ij 
are called the "(fc,m) spectral weight sum rules". The 
purely J = term in an "inclusive" sum rule (one with 
both J = and + 1 spectral contributions) is called 
"longitudinal" in what follows. 
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V us is extracted using flavor-breaking combinations 

SR w (s Q ) = [RZ(so)/\V ud \ 2 ] - [RZ (so)/\V us \ 2 } , (3) 

whose OPE representation, SRq PE (so), begins at dimen- 
sion D = 2, with a term proportional to m 2 . Experimen- 
tal determinations of Rudus( s o) over a range of Sq and 
w(s) then allow V us (and in principle also m s ) to be ex- 
tracted HHI1]. Explicitly, from Eq. ©, one has 



! RZ(so)/ 



R 



ud( S o) :r>w („ \ 

K^F " 6R ° pe(so) 



(4) 



The approach is favorable for the determination of V us 
because, at scales ^2 — 3 GeV 2 , and for weights 
used in the literature, SRq PE (so) is much smaller than 
the individual flavor ud and us spectral integrals (typ- 
ically at the few-to-several-% level). An uncertainty 
A {5Rope( s o)) in SRq PE (so) thus produces a fractional 
uncertainty (~ A (8R% PE {s )) /2R™ d (s )) in |K*| much 
smaller than that on SRq PE (sq) itself. Moderate pre- 
cision for SRq PE (sq) thus suffices for a high precision 
determination of \V US \, provided experimental errors can 
be brought under control [2j. At present, the accuracy of 
the method is limited by the sizeable (statistics limited) 
uncertainties on the us spectral distribution (which pro- 
duce ~ 3 — 4% uncertainties on the us spectral integrals, 
and hence 1.5 — 2% errors on | | ) . These errors will be 
radically reduced in the near future through analyses, al- 
ready in progress, of the new r decay data from BABAR 
and BELLE. We thus concentrate here on the issue of 
controlling OPE-induced errors. 



II. SOME TECHNICAL COMPLICATIONS 

The results reported below are all based on the ALEPH 
ud and us spectral data Q, ALEPH being the only col- 
laboration which has made its us covariance matrix pub- 
licly available. The reported results have been rescaled 
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to take into account current values of the various branch- 
ing fractions @ . In the case of the us data, this is done 
following the prescription described in Refs. 0. Details 
of the OPE input employed, and the method used to es- 
timate the D = 2 truncation uncertainty (which is much 
more conservative than those used previously in the lit- 
erature) may be found in Ref. [1]. 

The first important technical complication is the very 
bad behavior of the integrated longitudinal D — 2 OPE 
series. Not only does the series show no sign of converg- 
ing, even at the maximum scale, sq = m 2 , allowed by 
r decay kinematics @, Q but, even worse, it badly vi- 
olates spectral positivity constraints for all truncation 
schemes used in the literature (loj . Fortunately, the 
problem is easily handled phcnomcnologically. The ba- 
sic reason is that the longitudinal decay contributions 
are, for a combination of chiral and kinematic reasons, 
strongly dominated by the n and K pole terms [T3, El • 
These terms are very accurately known. The resid- 
ual non-pole ("continuum") contributions can, moreover, 
be well-constrained using the approaches described in 
Refs. d EH EI EI E3- Bin-by-bin subtractions of the 
longitudinal contributions to the experimental decay dis- 
tributions can thus be made, allowing direct determina- 
tions of the (0 + 1) spectral functions, from which FESR's 
not afflicted by the problematic longitudinal D = 2 OPE 
behavior can be constructed. The extremely bad behav- 
ior of the longitudinal OPE series makes inclusive anal- 
yses employing the longitudinal OPE representation un- 
tenable, and the longitudinal subtraction process unavoid- 
able. 

The second technical complication is the slow conver- 
gence of the £) = 2 J = + 1 OPE series, which is known 
exactly to O(a^) and given by [1] 
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where a = a s (Q 2 )/ir and fh s — m s (Q 2 ), with a s (Q 2 ) and 
m s (Q 2 ) the running coupling and strange quark mass in 
the MS scheme. Since independent high-scale determi- 
nations of a s (Mz) [a] correspond to a(m 2 ) ~ 0.10 — 0.11, 
the convergence at the spacelike point on |s| = sq, for 
s o < m r> 1S marginal at best. For the J = + 1, (k, 0) 
spectral weights, w( k ^(y) = (1 + 2y)(l - y) k+2 , with 
y = s/s = e**, the factor |1 - y\ k+2 = |2 sin{cj) / 2)\ k+2 
peaks more and more strongly in the spacelike direction 
with increasing k. Slow convergence of the integrated 
D = 2 series, deteriorating with increasing k, is thus 
expected for the (k, 0) spectral weights, an expectation 
borne out by the results of Refs. Fig. 1 of Ref. Q 

also shows no consistent common m s , V us fit region exists 
for the 0(a 4 )-truncated version of those (fc, 0) spectral 
weight sum rules employed previously in the literature. 

Since terms higher order in a s are relatively more im- 
portant at lower scales, a signal of the premature trunca- 
tion of a slowly converging OPE series is an unphysical 
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FIG. 1: so-stability plots for the (0,0) spectral weight, 
horizontal axis shows so in GeV 2 . 



The 



so dependence to the truncated sum. Such an unphysical 
dependence can be exposed by comparing the truncated 
OPE and corresponding experimental spectral integrals 
over a range of sq. This comparison is shown for the (0, 0) 
spectral weight, for a range of input m s = m s (2 GeV), in 
Fig. 1. The value of \V US \ needed for the spectral integral 
differences is obtained by fitting at s — m 2 . It is clear 
that no m s exists giving a good match between the OPE 
and spectral integrals. One should bear in mind that 
very strong correlations exist amongst the OPE integrals 
for different so and, similarly, amongst the spectral in- 
tegrals for different so- The results thus clearly point to 
problems with controlling OPE uncertainties in the (k, 0) 
spectral weight sum rules. That there is no obvious way 
to reduce the problematic D — 2 OPE truncation com- 
ponent of this uncertainty represents an important limi- 
tation to the use of the (0, 0) spectral weight sum rule as 
a means of determining | V us \ . 

The fact that |a s (Q 2 )| decreases as Q 2 = — s moves 
away from the spacelike point along the circle |s| = sq 
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allows one to improve the convergence of the integrated 
J = + 1, D = 2 OPE series by working with weights 
which emphasize those parts of the contour away from the 
spacelike point. Three weights, Wxo, Ww, and W20, de- 
signed specifically to produce improved convergence (in 
addition to improved control of other features of the OPE 
analysis) were introduced in Ref. [Til . A fourth such 
weight, w$, was introduced in Ref. [4J. The improved 
convergence, and resulting much improved sn-stability of 
\V US \, have been demonstrated previously 0, [j, [TT[ . Here, 
in Fig. [51 we show the analogues of Fig. 1 for the weights 
u>2o and w\o (the results for wiq and ws are intermediate 
in quality, and hence not shown explicitly) . The contrast 
to the (0,0) spectral weight case is immediately evident. 
We have also shown elsewhere (see Fig. 4 of Ref. [ij]) 
that an excellent common fit region for m s and | V us | ex- 
ists for the four non-spectral weights considered here, in 
sharp contrast to the lack of consistency found for the 
(k, 0) spectral weights. 

III. RESULTS AND PROSPECTS 

Single non-spectral weight, sq — rn^, fits, using as in- 
put the average, m s (2 GeV) = 94 ± 6 MeV, of strange 
scalar and pseudoscalar sum rule [T3| and Nf = 2 + 1 lat- 
tice results 15], advocated in the last of Refs. [l3|, yield 
\V US \ = 0.2209 ± 0.0029 e;C p ± 0.0017 t/l for w 20 , 0.22 1 ± 
0.0030 ea; p±0.0010 tfl for w w , 0.2206±0.0032 CTp ±0.0007 t/l 
for w w , and 0.2218 ± 0.0037 exp ± 0.0009 th for w$. The 
errors are dominated by the large uncertainties on the 
experimental us distribution. 

In view of the intrinsically slow convergence of the 
J = + 1, D = 2 OPE series, it is worth testing further 
the nominally improved convergence of the non-spectral 
weight sum rules. We do so by ignoring external infor- 
mation on m s and performing a combined, sq = rri^ 
fit for m s and \V US \, obtaining m s (2 GeV) = 96 ± 31 
MeV and \V US \ = 0.2208 ± 0.0052. The result for m s 
is seen to be in excellent agreement with the average 
of independent sum rule and unqucnchcd lattice deter- 
minations, already quoted above. This agreement justi- 
fies a combined sq — rr? T fit with the external average, 
m s (2 GeV) = 94 ± 6 MeV, as input. The result, 

\V US | =0.2209 ±0.0031 , (6) 

is compatible, within mutual errors, with both the re- 
cent Kta determination, \V US \ = 0.2249 ± 0.0019 
and the recent r^^/T^^] determination, \V US \ — 
0.2223^ ° 1 ° 3 26 [H. Eq. © represents our best determi- 
nation of \V US \ from present data. 

As for future prospects, re-running the combined non- 
spectral weight analysis with an assumed uncertainty of 
±5 MeV on m s (2 GeV) and us spectral errors reduced by 
a factor of 5 (a reduction expected to be readily achiev- 



able from the B-factory analyses), produces a total com- 
bined error < 0.0010 on \V US \. Further improvements to 
the theoretical component of this error are almost cer- 
tainly possible, once the currently rather large us errors 
above the K* are reduced, making it practical to con- 
struct, test, and work with additional weights having less 
strong suppressions of contributions from this part of the 
spectrum than do those employed above. 

We conclude by discussing two aspects of the cur- 
rent data situation which make it likely that higher cen- 
tral \V US \ values will be obtained from future analyses. 
The first point concerns the us data. At present, the 
branching fractions for observed strange decay modes go 
down only to the ~ 3 x 10~ 4 level (~ 1% of the total 
strange branching fraction) Q. Adding the contribu- 
tion of any additional, as-yet-undetected strange mode 
will necessarily increase \ V US \. The increase from a mode 
with branching fraction ~ 1 x 10~ 4 , for example, would 
be ~ 0.0003 - 0.0004. With the B-factory experiments 
already reporting preliminary results for strange mode 
branching fractions at the few-to-several 10~ 5 level [T3], 
it is clear that such missing mode contributions will no 
longer be a problem once the BABAR and BELLE anal- 
yses are completed. The second point concerns the ud 
data. With quoted errors, current experimental ud un- 
certainties produce an error of ~ 1/4% on \V US \. One 
should, however, bear in mind that significant discrepan- 
cies exist between the version of Py^J (s) measured in 
r decay and that implied by CVC and electroproduction 
(EM) data, even after all known isospin- breaking correc- 
tions have been taken into account [18j. The discrepan- 
cies are much larger than the nominally quoted errors on 
the r decay data, the direct r measurements and EM ex- 
pectations for the branching fractions of the Tr~ir°is T and 
7r~ ir~ 7r + ir° 'v T modes differing by 4.5er and 3.6a, respec- 
tively. Replacing the r tttt data with the nominally equiv- 
alent EM results would raise \V US \ by ~ 0.0017 - 0.0018. 
A similar replacement of the tt~ tt~ tt + tt° data would raise 
it by a further ~ 0.0016. While at present the r de- 
cay data is in better agreement than the EM data with 
certain OPE constraints [T{|, there is definitely room in 
the r decay data for somewhat lower ud branching frac- 
tions [20j . which would lead to larger \V US \. Resolving 
the t versus EM discrepancy is thus important for the 
determination of \V US \ using hadronic r-decay-based sum 
rules, as well as for resolving the status of the hadronic 
contribution to (g — 2) M in the Standard Model. 



Acknowledgments 

The financial support of the Natural Sciences and En- 
gineering Research Council of Canada is gratefully ac- 
knowledged. 



4 



OPE vs Spectral Integral (w ) 
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FIG. 2: so-stability plots for W20 (left panel) and W10 (right panel). The horizontal axis shows so in GeV 2 
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